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ABSTRACT 

A massive black hole binary might resonantly trap a star (e.g. a white dwarf) and the gas re- 
leased by its tidal disruption might emit electromagnetic wave signals around the coalescence 
of the binary. With post-Newtonian equations of motion including gravitational radiation reac- 
tion, we numerically studied resonant trappings by black hole binaries with mass ratio 1/100. 
It is found that 2:1 (and simultaneously 4:2) mean motion resonances of the binaries would 
be strong and could, in principle, draw small third objects deep into relativistic regimes {e.g. 
^ 10 Schwarzschild radii). The inclinations of the trapped objects could increase significantly 
and, in some cases, retrograde orbits could be realized eventually. 
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1 INTRODUCTION 

Various kinds of mean motion resonances are identified among 
satellites of Jupiter and Saturn (Murray & Darmott 2000). Their 
origins are explained by orbital evolution of the satellites due to 
their tidal interaction with the planets (Goldreich 1965). As for the 
orbital motions around the Sun, it is well known that Pluto and 
Neptune are in a 2:3 resonance (Murray & Darmott 2000). Mean 
motion resonances are also found in extrasolar planetary systems 
(Wright et al. 2011), including a large number of candidates re- 
cently provided by the Kepler mission (Lissauer et al. 201 1). Here, 
orbital migrations in protoplanetary disks would be the underly- 
ing process to yield the resonances (Snellgrove et al. 2001; Lee 
& Peale 2002, see also Raymond et al. 2008). Mean motion reso- 
nances, once established, are considered to play fundamental roles 
on evolution of multiple planet and satellite systems (see e.g. Mur- 
ray & Darmott 2000). 

Orbits of black hole binaries are shrunk by gravitational radi- 
ation reaction, and strong gravitational waves from merging mas- 
sive black hole binaries are promising targets of the Laser In- 
terferometer Space Antenna (LISA, Bender et al. 1998). One of 
the related phenomena in active debates is potential electromag- 
netic wave counterparts associated with the mergers. Seto & Muto 
(2010) discussed possibilities that stars (e.g. white dwarfs) might 
be resonantly trapped by massive black hole binaries and electro- 
magnetic wave signals might be produced by the gas released by 
tidal disruptions of the trapped objects (see Stone & Loeb 2010 
for similar arguments and also Fujii et al. 2010; Chen et al. 2011 
for possible resonant effects not associated with coalescences of 
binary black holes). In the paper they also examined the coorbital 
(1:1) resonances numerically and demonstrated that, depending on 
initial conditions, a test particle can be stably trapped by a black 



hole binary around its Lagrangian points L4 or L5 until the epoch 
close to the coalescence of the parent binary (see also Schnittman 
2010; Muto & Seto 2011). 

Stabilities of the coorbital resonances are relatively easy to 
study both analytically and numerically, but the critical aspect is 
how to originally settle the trapped objects around L4 or L5 un- 
der the presence of strong gravitational perturbations of the parent 
black hole binaries. In this paper, we analyze the standard mean 
motion resonances with ratios of the orbital periods m : n (m 7^ n) 
other than the atypical coorbital ones (with the ratio 1:1). While 
there are many parameters to set up the initial conditions for a 
triple system composed by a black hole binary and an additional 
small object, we here concentrate on the primary question whether 
a standard mean motion resonance can, in principle, trap a small 
object stably until a relativistic stage close to the merger of the par- 
ent binary. Below, by numerically integrating post-Newtonian (PN) 
equations of motion for triple systems, we will show preferable re- 
sults for 2: 1 (and simultaneously 4:2) mean motion resonances. 



2 DYNAMICS OF TRIPLE SYSTEMS 
2.1 Mass Ratios 

We study relativistic evolution of triple systems with masses 
mo, mi and m2, including only the gravitational interaction be- 
tween them. The particle mo(S> mi + 1712) is the central object 
and two point masses mi (inner one) and m2 (outer one) are orbit- 
ing around ?no. Hereafter we denote m = mo dropping the sub- 
script for simplicity, and put m = 1 in the geometrical unit with 
G = c = L In this subsection, we discuss how to choose the mass 
ratio m : mi : m2 for our numerical simulations. 

Using arguments of Hill stability for two particles mi and m2 
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orbiting around m in Newtonian dynamics, it can be analytically 
estimated how closely we can configure two (coplanar and nearly 
circular) stable orbits for a given mass ratio (mi + 7712 )/m (Glad- 
man 1993; Lee et al. 2009). Here we take (mi + m2)/m = 10"^ 
roughly corresponding to the stable limit for orbital periods of 2:1. 
For a black hole binary with a larger mass ratio (i.e. closer to 1), 
we need to increase the contrast of the two orbital periods (see also 
Holman & Wiegert 1999 and references therein) and the involved 
resonance should be higher than first order. But note that, once a 
mean motion resonance is realized, it can maintain stability for a 
higher mass ratio than the simple analytical predictions. For ex- 
ample, the allowed mass ratios (mi + m2)/m for coorbital (1:1) 
resonances become according to Hill stability, but a linear stabil- 
ity analysis around the Lagrangian points L4 and L5 provides the 
limit (mi + m2)/m < 0.04 (Gascheau 1843). 

For the internal weight between mi and m2, we examine the 
cases with mi < m2. This is because the orbital decay due to 
GW emission is faster for larger mass (see e.g. Landau & Lifshitz 
1971) and configuration with relatively approaching orbits is a key 
element for yielding capture into a resonance (Sinclair 1972; Hen- 
rard & Lamaitre 1983; Peale 1986). Hereafter we fix mi = 10"^ 
and m2 — 10^^, considering a fiducial astrophysical triple system 
composed by two BHs (lO^M© -f 10^ A/©) and a trapped IM© 
white dwarf. The numerical results in the next section are presented 
for these mass ratios, but we obtained qualitatively similar results 
for the combination IQ*^ Mq + 10* Af© + IMq. 



term H2PN is composed by many elements, and we use its lengthy 
expression given in Lousto and Nakano (2008) (see Appendix A in 
their paper). The 2.5PN term H2.5PN can be found e.g. in eq.(41) 
in Jaranowski & Schafer (1997). 

In this paper, we introduce new variables Sic, = Pia/nii to 
improve accuracy at numerically handling triple systems with large 
mass ratios (Seto & Muto 2010; Muto & Seto 201 1). This prescrip- 
tion is based on the technical reason, and we fully use the symmet- 
rical three body Hamiltonian without introducing approximations 
associated with ms <^ 1 (except for the numerical experiments 
discussed in Appendix A). In other words, the three particles are 
dealt equivalently in the post-Newtonian framework. Our equations 
of motion are given by appropriate partial derivatives as 

dxia _ 1 dH dsia _ 1 dH 
dt nii dsia ' dt rrii dxia 

The right-hand side of these equations are well behaved even in 
the test particle limit m^ — 0. The IPN contributions to the equa- 
tions of motion can be derived from eq.(3), but their explicit forms 
are also found in Lousto and Nakano (2008) (see their eqs.(4) and 
(5)). For the partial derivatives of the 2PN Hamiltonian, we gener- 
ate their fortran forms using Mathematica. The dissipative 2.5PN 
contributions for eq.(4) are presented e.g. in Galaviz & Bruegmann 
(2010) (their eqs.(12) and (13)). For our triple systems, the total 
number of equations becomes 18. They are integrated numerically 
as described in the next section. 



2.2 Equations of Motion 

To include general relativistic effects for triple systems, we use the 
2.5PN ADM Hamiltonian formally expressed as follows (Schafer 
1987; Jaranowski & Schafer 1997; Galaviz & Bruegmann 2010); 
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(1) 



where Hn, Hipm and H2PN are the Newtonian, IPN and 2PN 
tenns. These terms are conservative. H2.5PN is the 2.5PN dissipa- 
tive term due to gravitational radiation. In eq.(l), we introduced 
the set of parameters BpN = (6i,&2,fo2.5) for convenience to 
compare relativistic effects, and our main results are obtained for 
BpN = (1, 1, 1). In this paper, we negelct spins of the black holes. 

The ADM Hamiltonian H is originally given for the position 
variables xta and their conjugate momentums pia (i: label for the 
particles and a: for spatial directions). For example, the Newtonian 
term is written as 
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where we defined Pi = p^ ■ p^, dij = \xi 
(xi — Xj)/dij. The second and third summations in eq.(3) repre- 
sent interactions between different particles. Meanwhile, the 2PN 



2.3 Orbital Parameters 

We denote the coordinate distances di = |xi— xo|,(i2 = |x2— xo] 
and di2 = |xi — X2I between three particles mi,m2 and m(= 
mo). The semimajor axes ai and the eccentricities of the in- 
ner particle (i — 1) and the outer particle (i — 2) are numeri- 
cally determined from the consecutive maximum (apocenter dis- 
tance (1 + Bi)) and the minimum (pericenter distance (1 — e^)) 
of the radial separation di between the central object m and the 
particle m^. The orbital angular parameters are defined in a stan- 
dard manner as shown in Fig.l (see also Murray & Dermott 2000). 
The basic variables for studying resonant states are the longitude 
of ascending node Qi, the argument of pericenter uji and the mean 
anomaly Mi — 2n{t — Ti)/Ti {Ti. time of the pericenter passage, 
Ti: the interval between the passages). In the case of Newtonian 
elliptical motion. Mi is related to the true anomaly fi (shown in 
Fig.l) as 



= fi - 2e, sin fi + ^ef sin 2fi 



(3) 



0(e?). 



(5) 



3 NUMERICAL ANALYSES 

3.1 Initial and Termination Conditions 

We set up the initial conditions of our numerical calculations in 
the following manner. The outer particle m2(= 10^^) is placed at 
d2 = 300 with 62 — 0. The inner particle mi(= 10"^) is released 
around the pericenter at the distance di = SOOg^''^ with a control 
parameter q (q £ [0.27, 0.99] at 0.01 interval). We set its initial 
eccentricity ei ~ 0.15 and the relative inclination Ii ~ 0.01 [rad]. 
Although we do not go into details about the earlier histories of 
our triple systems, we expect that such an orbital alignment would 
be realized e.g. through interaction with the disk around the central 
black hole m. But, as we see later, exact coplanarity is not a critical 
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Figure 1. Definitions of orbital angular parameters. We take the reference 
plane (normal to the Z-axis) as the initial orbital plane of the outer par- 
ticle m2. The reference direction on the plane is the Jf-axis. The (green 
colored) objects are on the reference plane, and the (red colored) ones are 
on the (precessing) orbital plane of the inner particle mi . The longitude of 
the pericenter is defined by roi = Qi + uJi with the longitude of ascend- 
ing node Sli and the argument of pericenter lui shown above. The mean 
longitude is written by Ai = Mi + zui with Mi = /i + 0{ei). The 
parameters ro2 and Q2 for the outer particle m2 are defined similarly (for 
A2 see §3.3). During the evolution of the system, the particle m2 nearly 
stays on the reference plane with I2 ~ 0. 



requirement for yielding a stable resonant trapping. From Kepler's 
law and the relation di = ai{l — ei) for the pericenter distance, the 
ratio of the initial orbital periods Ti and T2 is given as Ti /T2 
(ai/a2)^''^ - 1.2q. 

We define the reference plane for the orbital angles by the 
initial orbital plane of the outer particle 7712 (see Fig. I). Given 
nil <^ m2, the total initial angular momentum of the triple is 
nearly normal to the reference plane (parallel to the Z-axis in 
Fig.l). The two particles mi and m2 are initially set nearly at con- 
junction with the argument of pericenter cji ~ 7r/2. 

Next we describe when we stop our runs. Given our pertur- 
bative (post-Newtonian) treatment of the relativistic effects, we 
conservatively terminate the integration when one of the follow- 
ing close encounters happens for the first time; (I) di < 10m, 
(II) di2 < IO7TI2 and (III) d2 < 10m. We set another termi- 
nation condition (IV) di > 30d2, as a simple ejection of the 
lightest particle mi from the system. These criterions are se- 
lected primitively for particle systems. But the tidal radius of 
a white dwarf (mass IA/q and radius 10*km) is given by ~ 
lO''(m2/lO^M0)^/^km for a black hole of mass m2, and larger 
than 10m2 ~ lO''(m2/lO^Af0)km used for the condition II. 
Therefore, the white dwarf would be tidally disrupted before the 
condition II is applied. On the other hand, the central black hole m 
can directly swallow a white dwarf. 

We solve the equations of motion (4) using a fifth order 
Runge-Kutta scheme with an adaptive step size control (Press et 
al. 1996). Here we briefly describe two of the numerical tests done 
for our code. These are performed for the specific sets of parame- 
ters (a2, 62) = (100, 0.1) and (50,0.1). 

(i) Binary evolution by gravitational radiation reaction. For a 
purely binary (tentatively putting mi =0), we compare our numer- 
ical results da2/dt and de2/dt with the corresponding analytical 
expressions obtained by Peters (1964). Since the latters are given 



for BpN = (0, 0, 1) in our post-Newtonian parameterization, we 
prepared the numerical results in the same setting Bp jv — (0,0,1). 
We found that the differences between the analytical and numerical 
ones are less than 1%. 

(ii) Conservation of Hamiltonian for triple systems. For sev- 
eral resonant triple systems with mi 7^ and Spjv ~ (1, 1,0) 
(namely excluding the dissipative 2.5PN terms), we checked con- 
servation of the total Hamiltonian H = Hn + Hipm + H2pn- 
The variation rate dH/dt is less than 10~^ of the expected energy 
loss rate due to the 2.5PN terms. 



3.2 Final States 

Our primary interest in this paper is stable resonant evolution of rel- 
ativistic triple systems. Therefore, we first provide the final semi- 
major axes a2,fin of the outer particles m2 at the end of our cal- 
culations (started from a2Ani = 300). In Fig.2, the results a2.fin 
are shown with the filled symbols. These are obtained with the pa- 
rameters BpM ~ (1,1,1). The squares are for the case I (close 
encounter between the inner one mi and the central object m), the 
circles for the case II (close encounter between mi and m2), and 
triangles for the case IV (ejection of mi). No run ended with the 
condition III. Fig.2 shows that, for q — 0.33 ~ 0.38, the inner 
particles mi were stably trapped by the outer ones m2 down to rel- 
ativistic regime 02, /in ~ 20. These are due to 2:1 mean motion 
resonances. In the next section, we pick up one of them (q — 0.33) 
and follow its evolution. Another dip around q = 0.55 was caused 
by 3:2 resonances that exist from the beginnings of the runs (as 
expected from Hill stability). 

In Fig.2 we also show the inclination sinJi around ~ 500 
cycles before the end of our calculations. The crosses are for retro- 
grade orbits with |Ji| > 7r/2 and the pluses are for prograde ones 
with |/i I < 7r/2. While we always had /2 ^ for the outer parti- 
cle, the angle Ii could grow significantly. In particular, for deeply 
trapped ones (with small 02, /in), the inner orbits could even be- 
come retrograde | Ji | > n/2. Many of them ended with the condi- 
tion I. But the condition II is realized in some cases (e.g. q = 0.40 
with a2jin ~ 60), and, for the fiducial system containing a white 
dwarf, we might observe electromagnetic wave signals associated 
with the tidal disruptions. One the other hand, in our numerical 
samples, the ejections of mi were common outcomes for the 3:2 
resonances. 

In addition to the standard setting BpN = (1,1,1) men- 
tioned so far, we examined evolution of the triple systems with 
BpN = (1, 0, 1) and (0, 0, 1) to study relativistic effects. In this 
comparison, the three particles are again handled equivalently with- 
out using approximations based on mi <^ 1, and the dissipative 
2.5PN term was remained to generate orbital migrations. The re- 
sults are presented in Figs. 3 and 4. We found that the overall be- 
havior of Fig.2 is unchanged for the IPN case given in Fig.3 with 
BpN = (1,0, 1). The minimum value of the final semimajor axis 
0,2. fin ~ 20 is similar to the 2PN results, and retrograde orbits 
are also observed. In contrast, the Newtonian results in Fig. 4 with 
BpN = (0,0, 1) are quite different from other two. For exam- 
ple, we no longer have retrograde inner orbits. While the 2: 1 res- 
onances are realized at lower q, compared with the 2PN and IPN 
cases, the minimum value of the final semimajor axes increases to 
0,2. fin ''^ 75 in Fig.4. These results indicate that the relativistic 
effects help the long term resonant trapping of binary black holes. 
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Figure 2. The final states of the triple system with masses m = 1 
(central object), mi = 10~^ (inner particle) and m2 = 10~^ (outer 
particle). The results are obtained with the full 2.5PN equations of mo- 
tion [(fci, fe, fe.s) = (1,1,1)]- The initial semimajor axis of 7x12 is 
o.2,ini = 300 and the ratio of initial two orbital periods is T\/T2 ~ l.2q. 
The filled symbols represent the final semimajor axes 02, /in [squares for 
di 10m (case I), circles for di2 ^ lOmi (case II) and triangles for 
ejection with di > 30(i2 (case IV)]. The crosses and pluses on the solid 
curve show the inclination sin /i at ~ 500 rotation cycles before the ter- 
mination of our integrations (crosses: retrograde orbits with > 7r/2, 
pluses: prograde with |/i| < 7r/2). 





1.0 




0.8 








0.6 


'S 




a 


0.4 


c 

<^ 




a 


0.2 




0.0 







( ' 'VA 
• 


A«AAA«AA 


' AOTAAWA«AAAA«AAAi 










A*** 






• 




A 






• 








"M. ■ ■ 


IT 










H 1 l-Ll 


1 1 1 1 1 1 1 1 


, HH 1 1 ill 1 lU li 1 li 1 il 1 1 1 li 



0.3 



0.4 



0.5 

q 



0.6 



0.7 



Figure 3. Same as Fig. 2 but the 2PN term is excluded [namely with 
(fel, 62, fe.s) = (1, 0, 1)]. The overall profiles of the semimajor axes and 
the inclinations are similar to the results in Fig. 2 obtained with 2PN terms. 



3.3 Evolution of Orbital Parameters 

Here we pick up the full 2.5PN mn for the initial parameter 
q = 0.33, and examine its time evolution. Our integration ended 
at a2,fin = 22 where the retrograde inner particle mi has a close 
encounter with the central one (the termination condition I). This 
run is a typical example for the long-term 2: 1 resonances appeared 
in Fig.2 around 5 = 0.35 with the filled squares. For reference, we 
also show the corresponding Newtonian result (with q — 0.33) in 
Fig.6 

Since the orbital decay is determined by gravitational radiation 
and we have the simple relation m ^ m2 S> mi, the semimajor 
axis a2 is related to the remaining time tr and the rotation cycles 
Nr before the coalescence of m and m2 as (Peters 1964) 

256 m2 V m / 647rm2 \ m / 



Figure 4. Same as Fig.2 but the IPN and 2PN terms are excluded [namely 
with (fei, ti2, fe.s) = (0,0,1)]. We no longer have retrograde orbits 
(crosses appeai'ed in Figs. 2 and 3) and minimum value of the final semi- 
major axes increased to a2,/i„ ~ 75. 



(ort^ ~ 3.1(a2/100m)''yrandiV,. ~ 5.0 x 10''(a2/100m)5/2 for 
the fiducial triple with m = W^'AIq and m2 = lO''i\/0). There- 
fore we use the semimajor axis 02 of the outer particle to describe 
the evolutionary stages of the system. The basic results are summa- 
rized in Figs. 5 and 6 with time elapsing leftward from 02 — 300. 
The PN effects largely change the evolution of the eccentricity and 
the inclination of the inner particle. Below, we mainly comment on 
the PN calculations given in Fig.5. 

We show the semimajor axis ai of the inner orbit in Fig.5a. It 
initially oscillates until the resonant capture around 02 ~ 285 (as 
confirmed later through a resonant argument), and then smoothly 
shrinks as ai ^ 2~^''^a2. In Figs. 5b and 5c, we show the eccen- 
tricities ei and 62 of the two particles. They show short period os- 
cillations. The inner one becomes ei ~ 0.2 around 02 — 250, 
while the outer one keeps small values 62 < 2 x 10^^. 

The relative angle SI2 — Sli is presented in Fig.Sd. It soon 
settles to Q2 — ^1 ^ 71". This means that the components of the an- 
gular momentums of mi and m2 projected to the reference plane 
are nearly in the opposite directions. The oscillation in the early 
stage (with small I/2I), would be partly due to our choice of the 
reference plane. For the outer particle, we initially have J2 = 
and the longitude of the ascending node f72 is not well defined. But 
other angles Ai , A2 , f2i , zui and VJ2 can be determined almost un- 
affected by this choice, as understood from their intrinsic geomet- 
rical meanings. Since we have I/2I <C 1 and 62 <C 1, we determine 
the mean longitude A2 using the time when m2 pass the XZ-plane 
in Fig.l. With mi ^ m2, the initial total angular momentum is 
nearly normal to the reference plane. 

In Fig.5e, we show the inclination angle Ji. The magnitude 
of the angle I2 for the outer one stays at I2 ~ 0, as expected 
from mi <C m2. At the Newtonian order, the angular momen- 
tum Ji of the particle m^ is related to the orbital period Ti as 
r'''^l-e?)i/2. Therefore, we define the combination 



E = y^l - ej sinJi - 1000 ■ 2^''^^1 - e2sinJ2 



(7) 



for the component of the total angular momentum projected to the 
reference plane, and plot it in Fig.5e. The result E ~ together 
with the relation ^2 — f^i ~ tt shows that the orientation of the 
total angular momentum is nearly parallel to the Z-axis and almost 
unchanged. A similar behavior is observed in the Newtonian results 
(see Fig.6e) and also in numerical simulations of resonant planet 
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Figure 5. Evolution of the triple system with q = 0.33 and Spjv = 
(1, 1, 1). The horizontal axis is the semimajor axis a2 of the outer particle 
from the initial value a2^ini = 300 down to the final value a2,/in = 22. 
The calculation is terminated by the condition II (close encounter between 
mi and m2). The panel a) shows the semimajor axis of the inner particle, b) 
and c) The eccentricities of the inner and outer particles respectively, d) 
The relative angle f!2 — f2i between the longitudes of ascending nodes, e) 
The inclination sin Ii and the combination S defined in Eq.JT) for the total 
angular momentum. The inner orbit becomes retrograde at a2 < 60. f), g) 
The resonant arguments 8ei = 2A2— Ai— roi andS/i/2 = 2A2— Ai— Hi 
plotted once every 500 rotations of m2 . h) The relative angle 1172 ~ zui 



migrations (Thommes & Lissauer 2003). These adiabatic results 
are reasonable, since the time scale of the precessions Qi is much 
shorter than that of the orbital decay. 

In Figs.5f and 5g, we present the resonant arguments 



2A2 — Al — Ci7l, 



4A2 - 2Ai - 2fii 



(8) 



defined for a 2:1 eccentricity resonance and a 4:2 inclination res- 
onance respectively. Due to the symmetry with respect to the ref- 
erence plane, the inclination resonances start from the second or- 
der in the standard classification of resonances (Murray & Darmott 
2000). As shown in panel Fig.5f, the argument 6ei initially cir- 
culates, but turns to librate around 6^1 ~ at 02 ~ 285 where 
the inner particle is captured into the 2:1 eccentricity resonance. 
The libration amplitude gradually decreases, and the inclination Ii 
shortly increases up to sin Ii ~ 0.2. On the other hand, throughout 
this run, the relative angle UJ2 — yji and thus the resonant argument 
6^2 = 2A2 — Al — 1^2 circulates (see Fig.5h). But for some other 
initial parameters q {e.g. g = 0.36 with the final values {a2ji„ , Ii ) 
similar to g = 0.33, see Fig.2), we found structured distribution 
for the combination uj2 — ^1. Our Newtonian njn in Fig.6h also 
shows the libration VJ2 — yJi ~ tt at 02 J; 280. 

At a2 ~ 265, in addition to 6^\, the argument On is cap- 
tured into libration around 9n/2 ~ 37r/2. From Figs.Sd, 5f and 
5g, we can see that other two arguments, 4A2 — 2Ai — 20,2 and 
4A2 — 2Ai — 0,1—0,2 are also in libration states, and we have 
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Figure 6. Same as Fig.5 but without the IPN and 2PN terms (Bpjv 
(0, 0, 1)). This run ended with the condition 1 around 02 ~ 80. 



\vD\ — fii I ~ 7r/2 as in the case of the Kozai libration (Kozai 1962) 
that is caused by secular terms and known to be vulnerable to rel- 
ativistic effects in the course of orbital migration. In Fig.5g, the 
distinct outlier point around 02 = 110 is associated with ei <c; 1 
for which the mean longitude Ai would not be well-behaved. In the 
Solar system, two Satumian satellites, Mimas and Tethys are in the 
inclination resonance with the argument 4A2 — 2Ai — Oi — O2. 

ConcuiTently with the libration of On, the inclination I\ starts 
to increase due to accumulation of resonant effect. Similar evolu- 
tion appeared in the 3:2 resonances. Growth of inclination accom- 
panied by simultaneous (eccentricity and inclination) resonances 
is observed also in numerical simulations of planet migrations in 
Newtonian dynamics (Thommes & Lissauer 2003). But, there, rel- 
atively large eccentricities ei ~ 0.6 are required to yield the si- 
multaneous resonances. Interestingly, our Newtonian simulation in 
Fig.6 shows that we indeed need large eccentricities ei ~ 0.6 for 
the onset of the simultaneous resonances around a2 ~ 160. 

Now we briefly discuss the differences between the PN and 
Newtonian calculations. Here, it is important to note that the si- 
multaneous resonances {6^\ and 9n, see eq.(8)) are realized with 
the relation 0\ ~ roi. First we examine the Newtonian run given 
in Fig.6. In the eccentricity resonance (with the argument 6ei) for 
a nearly coplanar Newtonian system with ei ^ 1, the inner peri- 
center retrogrades as > uji oc 1/ei (Murray & Dermott 2000). 
In Fig.6, the simultaneous resonances (&ei and On) are established 
at a2 ^ 160 when the eccentricity ei reaches 0.6. In the earlier 
stage 02 > 160, we found > fii > zui, reflecting the smaller ec- 
centricity ei and the resultant higher regression rate uji. Thus the 
condition Oi ~ rbi is first met at the relatively large eccentricity 
ei 0.6. 

Next we analyze the PN calculations given in Fig.5. The PN 
effects generate a positive contibution {vji) pn > to the apsidal 
precession (Landau & Lifshitz 1971). Meanwhile, from close com- 
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parisons between the runs in Figs. 5 and 6, we found that the nodal 
precession rate Ct\ is less affected by the relativistic effects for the 
nearly coplaner triples. Therefore, the difference zui — Q.\ < Q 
for the Newtonian results at ei < 0.6 is compensated by the PN 
correction (Ti7i)pjv > 0, and the condition -riji ~ Q,\ (thereby in- 
crease of the inclination I\) can be realized at a smaller ei, as seen 
in Fig. 5. 

During the stable resonant evolutions of the two cases, close 
encounters between mi and m2 are avoided. The inner eccentric- 
ity ei shows remarkable contrast between Fig. 5b and Fig.6b, and 
this is an important parameter for the condition I, as we see be- 
low. In Fig. 5, the inner particle is on a polar orbit relative to the 
outer one around a2 ~ 60, and then retrogrades with the final value 
sin/i ~ 0.5 (7i ~ 150°). Just before the end of this run, the ec- 
centricity ei goes up (not sufficiently captured in Fig. 5b), and the 
calculation is ended by the condition I with the pericenter distance 
(1 — ei)ai ~ 10. In contrast, for the Newtonian sytem in Fig. 6b, 
the inner eccentricity ei increases almost monotonically, and the 
condition I is met earlier at 02 ~ 75. 



4 DISCUSSIONS 

We have studied relativistic evolution of the resonant triple system 
with the post-Newtonian approach including terms up to 2.5PN or- 
der. We mostly fixed masses of the system at m = 1 (central ob- 
ject), mi = 10^^ (inner particle) and m2 — 10^^ (outer parti- 
cle). Here we briefly summarize our basic results. We found that 
an eccentricity-type 2:1 resonance (together with inclination-type 
4:2 resonances) is strong and the outer particle could, in principle, 
trap the inner one down to the relativistic regime with final distance 
02 ~ 20. Here, the PN effects play important roles for the resonant 
evolution. For example, they help the onset of the simultaneous 
resonances at a small eccentricity ei. For the PN system, the inner 
eccentricity ei can remain small, and the close encounter between 
m and mi is delayed, compared with the corresponding Newtonian 
simulation in which the eccentricity ei continues to grow. During 
the simultaneous resonances, we observed rising of the mutual in- 
clination angle, and even retrograde orbits could be realized. 

Next we describe astrophysical implications of our results, 
particularly for gravitational wave astronomy. As an example, we 
follow the evolution presented in Fig. 5 and introduce relevant phys- 
ical units by taking m — IQ^ Mq. The resonant trapping of the in- 
ner one mi continues down to a2 ~ 22, where the system becomes 
unstable and a close encounter between mi and m occurs with dis- 
tance d\ ~ 10m. At this point, the gravitational wave frequency 
from the main binary m-m2 is 2/T2 ~ 6mHz. The time before 
its coalescence is tr ~ 2.7 days and we have 2Nr = 2300 grav- 
itational wave cycles left. Before this transition, the inspiral wave 
of the main binary has a phase correction of C'(mi/m2) due to 
the trapped particle mi, compared with a purely two body system 
m-m2. It might be possible to effectively describe the waveform 
of the triple system in a simple manner with a small number of 
additional fitting parameters. Otherwise we might not make a long 
term coherent integration of the gravitational wave signal before 
the transition. 

If the trapped particle mi is a neutron star (or a stellar mass 
black hole) and plunges into the outer black hole (as for q = 0.40 
in Fig. 2), the associated gravitational wave signal might be detected 
by future ground-based interferometers as a precursor to the merger 
of the parent binary that might be observed with space interferom- 
eters. But these compact objects would also have chances to fly 



out from the systems at very high velocities > 0.1c. In the case 
of a white dwarf, we might observe characteristic flares of electro- 
magnetic waves due to the gas released by its tidal disruption (Rees 
1988). Detections of such electromagnetic wave signals would have 
significant impacts on astronomy and cosmology (e.g. identifica- 
tion of host galaxies). 

Here, we comment on possible expansion of our study. While 
we have mainly investigated triple systems already captured in first- 
order resonances, it would be meaningful to analyze the earlier evo- 
lutionary stages. The later stages around the binary coalescence 
should be also studied, including strong gravity beyond our per- 
turbative treatment. Spins of black holes might become important 
there, especially for highly inclined orbits. In addition, various as- 
trophysical processes {e.g. interaction with a gas disk, tidal defor- 
mation of extended stars) would be worth studied. 

This work was supported by ISPS grants 20740151 and 
22.2942. 
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APPENDIX A: COHERENT RADL4TION REACTION 
FORCE ON THE SMALL BODY 

When two particles mi and m2 move incoherently around m, 
the frictional forces caused by gravitational radiation can be re- 
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garded as independent. However, in the case of resonant motions, 
there can be coherent frictional effect between them. For an ec- 
centric inner particle in a 2: 1 resonance, we can expect this in the 
quadrupole order of gravitational radiation, as indicated by cou- 
pling between quadrupole moments of the two particles. Using the 
approximate formula for radiation reaction force given in Landau & 
Lifshitz (1971), we roughly estimated the frictional work imposed 
on the inner particle mi by the quadrupole moment of the outer one 
m2(S> mi), and found that this work alone would not be sufficient 
to keep the resonant state with respect to the infalling outer one. 

For reference, we also made the following numerical exper- 
iment in order to shut down the coherent frictional force caused 
by the main binary on the small inner particle. In eqs.(4), we put 
B2.5PN ~ for the lightest body (i = 1) but keep B2.5PN ~ 1 
for other two bodies (i = and 2). Then we obtained numerical 
results similar to Fig. 2. 

The standard resonant transfer of energy and angular momen- 
tum from mi to m2 would be the dominant and sufficient process 
to keep the stable resonant evolution as seen in Fig.2. In the case of 
a coplanar triple with £2 <^ 1, such transfer would be realized by 
conjunctions after the pericenter passages of the inner particle (see 
e.g. Murray & Dermott 2000). The pericenter are known to retro- 
grade in Newtonian dynamics, and its rate is proportional to oc e]~^ 
(Murray & Dermott 2000). 



